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Production of Bloch waves during a rapid quench 
is studied by analytical and numerical methods. The 
density of Bloch waves decays exponentially with the 
quench time. It also strongly depends on temper- 
ature. Very few textures are produced for temper- 
atures lower than a characteristic temperature pro- 
portional to the square of the magnetic field. 

Topological defects appear in many condensed mat- 
ter systems They are also believed to play a role 
in cosmology [|^ and in nuclear physics ||]. Kibble ^] 
pointed out that the defects can be produced in signifi- 
cant numbers during a rapid second order transition from 
the disordered to the symmetry broken phase. The speed 
of light was the key factor in his estimate of the defect 
density just after a quench. Zurek Q stressed the impor- 
tance of the nonequilibrium Ginzburg-Landau dynamics 
of the order parameter. According to the Kibble-Zurek 
scenario the density of kinks, vortices or monopoles scales 
like a certain defect-dependent power of the quench rate. 
The theory was verified in several experiments . 

Both theory and experiment concentrated so far on 
hedgehogs, which have a nontrivial asymptote at spatial 
infinity. Textures are a wide class of topological defects 
with a constant asymptote. A texture which arises in a 
D + 1 component spin system has D spatial dimensions. 
As the spin vector is constrained to have definite mag- 
nitude, it lays on a I? dimensional sphere. An external 
magnetic field prefers, say, downward spin orientation. 
A spin field has finite energy if its spins point down at 
spatial infinity. The D dimensional space is compactified 
to a D dimensional sphere. A finite energy spin field is 
a map — > S^. If the map has a nonzero winding 
number, then it can not be continuously deformed to a 
uniform spin-down ground state. 

The simplest example of a topological texture is a 
Bloch wave in one spatial dimension. Its spin is an in- 
plane arrow of definite length, which points down far from 
the wave. As we pass through the (anti-)texture from the 
left to the right the spin makes the full (anti-)clockwise 
rotation. The Bloch wave has a definite size as a result 
of balance between the Zeeman and the exchange energy. 

In this paper we study the production of Bloch waves 
during a rapid quench. We believe that our results 
are qualitatively correct for all topological textures like 
skyrmions in planar ferromagnetic systems (e.g. quan- 



tum Hall ferromagnets magnetic bubbles or Hopfions 
in three dimensional ferromagnets |^ and baryons in the 
Skyrme model Q . In addition, Bloch waves are formally 
equivalent to fluxons in long Josephson junctions 

The Bloch wave is topologically stable because its 
spin is constrained to lay on a circle. Above the criti- 
cal temperature the spin (the order parameter) expecta- 
tion value is zero. There are no textures there. Some 
nonzero density of textures can be created during a 
paramagnet— >ferromagnet transition. No textures are 
created in a quasistatic transition because the correlation 
length diverges at the critical point and the ferromagnetic 
phase is entered already with a long range order (quasi 
long range order in one spatial dimension). However, if 
the transition preceeds at a finite rate, then, according to 
the Kibble-Zurek scenario, the system passes the critical 
point out of thermal equilibrium. It enters the ferromag- 
netic phase with a finite correlation length. This lack of 
long range order can manifest itself with a finite density 
of topological textures. We will see that the density of 
Bloch waves created in this way depends not only on the 
transition rate but also on the external magnetic field. 
The leading dependence of the density on the transition 
rate is exponential, which is much in contrast to the power 
law dependence characteristic for hedgehogs. 



Time Dependent Ginzburg-Landau Model 

Let us define a one dimensional time dependent 
Ginzburg-Landau model by the stochastic field equations 
for the two spin components 4>a(t, x), a = 1,2, 



(1) 



where ' = dt, ' = dx- r]a{t,x)''s are Gaussian noises of 
temperature T with correlations 

{Va{t,x)) = , 

{r]a{ti,Xi)rif3{t2,X2)) =2TSal3S(ti~t2)S{xi~X2) ■ (2) 

The coefficient a{t) is time dependent. We consider a 
symmetric linear quench 



if t < 

if < i < 2r 

if 2T<t 



(3) 
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Before the quench, for i < 0, the system is in a symmetric 
(paramagnetic) phase (a > 0), during the quench, at 
t = T, it undergoes a transition from the symmetric phase 
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{a{t < r) > 0) to a broken symmetry (ferromagnetic) 
phase {a{t > r) < 0). FinaUy it settles down at a{t) = 
-1. 

In the final symmetry broken phase (a — —1) the spin 
field tends to live at the bottom of the Mexican hat po- 
tential V = -|(^| V2+|0|'*/4. The bottom of the potential 
is the circle |(/)| = 1. One can consider fluctuations in the 
direction normal to this valley. Their characteristic first 
relaxation time is 0(1). The potential plus the exchange 
energy generate a length scale 0(1) for these normal fluc- 
tuations. A nonzero magnetic field B removes degeneracy 
along the circle. We assume that B is small so that we 
can approximately write 



h + ii 



1)2 = e*x+^- 



(4) 



for the field at the bottom of the potential valley. The 
spin length is a hard mode so it has a fixed magnitude 
but its orientation (phase x), which is a soft mode, can 
vary in space and time. Eqs.(^ with a = — 1 and 77 = 
reduce to the effective equation 



X = X 



B sin X J 



(5) 



which can be recognized as a diffusive sine-Gordon equa- 
tion. The length scale in this effective model is B^^/^ 
and the second relaxation time is B~^ . If, as is generic, 
B « 1, then these two scales are much longer than 
those characteristic for the hard magnitude mode. The 
model (|5|) has a static soliton solution 

B^/'^x 

x(a;) = — 4arctan[tanh( — ^ — )] mod 2tt . (6) 
This soliton is just the Bloch wave. 



Instantaneous Quench 

Once we distinguished between hard and soft modes 
with corresponding length and time scales, we can set 
to the instantaneous quench (r = in Eq.(|3|)). In the 
symmetric phase (t < r) the order parameter fluctuates 
around the ground state (0i, ^2) = (0, —B). li B « 1 
and the temperature T is moderate, then the qubic non- 
linearity on the RHS. of Eq.(|l]) can be neglected. We are 
in the Gaussian regime. The system before the quench 
has a unit correlation length. As the quench is instan- 
tenous, it enters the symmetry broken phase with this 
correlation length unchanged. Its further evolution can 
be unambigously divided into two stages. 

The first stage is very short, it lasts for 0(1) units of 
time. This is the time the order parameter needs to roll 
down to the bottom of the sombrero potential. At the 
end of this stage the spin already has a fixed magnitude, 
101 « 1. The phase x is still chaotic but it is correlated 
over the length scale 0{1). The field already contains 
Bloch waves but they are very thin (0(1) width) and 
distorted, see the plot in Fig.l. 



In the second long stage, which lasts for another 
0{B~^) units of time, the initial "baby Bloch waves" 
grow in size, according to Eq. (^ , until they become fully 
pfledged Bloch waves (H), compare Fig.l. 

Wc want to estimate the density of the Bloch waves 
after the quench. Let us first estimate the density but of 
the baby Bloch waves. We define the position of the (anti- 
)baby Bloch wave by a point x where the spin points up. 
In the gaussian approximation the field fiuctuates around 
((/)i,02) = (0,-5) at i = 0~. The quadratic potential 
is turned upside down at t = 0. The field at x rolls 
down the slope of the sombrero potential to the upwards 
spin orientation iff 4'2{x) > B and (/'i(a;) = 0. is 
an average of (j) over the unit correlation length. In our 
calculations below we take this average by introducing 
an ultraviolet cut-off in momentum space at |fc| = 1. 
In the gaussian approximation the probability that the 
two conditions hold simultaneously is a product of their 
probabilities. Let us work out the probabilities one by 
one. 

The density of the points such that 0i (0, x) = can 
be worked out with the general formula |lffl 



N[(l)i = 0] = - 



1 {4>[{o,x)4>[{o,x)) 



(0i(O,a;)0i(O,x)) 



0.17 



(7) 



In the second equality the equilibrium two point corre- 
lation function in momentum space at t = 0~ was taken 
into account. 

Now we estimate the probability that 4'2{0~,x) > B. 
At t = 0" (f>2 fluctuates around —B. The magnitude of 
the fluctuation squared can be found as 



g^ = {[^2{Q-,x)+B]^) = - 



dk 



1 



1 



T 

F^2 



The normalized gaussian probability distribution of 

^2(0^, a;) isp(02) = exp-(02 + S)V2gVv/2^- With 
this distribution the desired probability is 



>B]^ 



+B 



) 1 9 jp 

p[x) = -{l-Erf[—]}, (9) 



where we use Mathematica's Erf[x] — exp(— 
The density of baby Bloch waves scales like 



2B. 



Nbaby = N[^i = 0] X P[02 >B]= 0.08 ci {1 - Erf[c2^]} 

(10) 

with extra parameters ci, C2 — 0(1) to reflect the order- 
of- magnitude nature of our calculations. 

Nbaby does uot distinguish between textures and anti- 
textures. Some pairs of baby textures and antitextures 
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will annihilate. In the second slow stage of the relaxation 
process a baby texture grows in width to its full mature 
size « 3 X If there is a baby antitexture within a 

comparable distance, the texture and antitexture anni- 
hilate each other. An average topological charge density 
within this distance is zero. Fluctuation around this av- 
erage density is the final density of Bloch waves N and 

it is proportional to \/ x A^; 



baby 



N = OmCi B^^^Jl-Erf[C2^] N « B^/^ 



(11) 

with constants Ci,C2 = 0(1). The formula is accurate 
in the dilute regime, N « B^^^. When extrapolated 
to the oversaturated regime, it gives an upper estimate. 
The curve ( pl| ) favourably compares with the results of 
our numerical simulations, see Fig2. 

Density of Bloch waves depends in a critical way on 
temperature. Very few textures are produced for tem- 
peratures lower than B^. This can be easily understood. 
If T << B^, then just before the quench the field is lo- 
calized in a very close neighbourhood of the ground state 
((/>!, 02) = (0, — -B). After the quadratic potential is re- 
versed, the field uniformly (up to negligible fluctuations) 
rolls down to (0, —1) without any chance to wind around 
the top of the Mexican hat potential. 



Finite Rate Quench 

We first estimate the density of baby Bloch waves 
at the end of the quench at t = It. Let us define 
rait) — {(/)2{t,x)) as a noise average of (j)2 or equivalently 
(j)2. A noise average of the linearized a = 2 component 
of EqsM is 



m{t) = -a{t) m{t) - B 



(12) 



The solution with the initial condition m(0) = —B and 
a{t) given by (||) at the end of the quench is given by 



m(2r) 



-B - bV2^ Erf[y/7/2] 



T»l 



(13) 



At the end of the quench the fluctuations of 02 around 
its average m(2T) are 



[ErfiVTik' + 1)) - Erf{V^{e ~ 1))]} 



T»l,k=^l 



(14) 



For r >> 1 the fluctuations tend to 51^ « T e'^/2\/2. 
Similarly as for the instantenous quench the probability 
that 02 (2t, x) > B is given by 



P[02(2t) > B] 



Ir r-B-m(2T),, 

-{l^Erf\ } 

2 x/2^ 



T 



8^/27^2rB2 



exp{-^^^^ } . (15) 



At the end of the quench the density of the points such 
that 0i(2t, x) = is given by 



N\6 =01 = - / (^i(2T,x)0'i(2r,x)) 
^ ttV (0i(2T,a;)0i(2T,x)) 



JZ^ dk fc2 G(2r, T, k) 
\\ dk G(2r, T, k) 



27rV^ 

The density of baby Bloch waves is 



(16) 



Nbaby{T»T/B^) 



T 



32V27r^TB^ 



, 2^/27^TB^ 
exp{ } 



T 



(17) 



We obtain the density of Bloch waves in a similar way as 
for the instantaneous quench, 



NiT»T/B')^]J^xNtaby = 

= 0.035 i?i(-|-)i/4exp{-i?2^^^^} , (18) 

with constants Di,D2 = 0(1). The density of textures 
decays exponentially with r on the time scale 0{T/B^). 
Eq.(|l8[) is compared with numerics in Fig. 3. 



Conclusion. 

Almost no textures are produced for T less than B^. 
Their density decays exponentially with the quench time. 
The time scale for this exponential decay is of the order 
of the first relaxation time. 

Our findings have qualitative implications for higher 
dimensional textures. In the case of Skyrme solitons 
which model baryons, an effective "magnetic field" is pro- 
vided by the pion mass term in the Lagrangian. The 
density of baryons just after a chiral transition should 
be supressed as compared to the density predicted by 
the standard Kibble-Zurek scenario. This case deserves 
a separate quantitative study; the study is given in a 
preprint UM- 
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This work extends the Kibble-Zurek scenario to the 
case with a bias field which exphcitly breaks the 0{D + 1) 
symmetry. In the standard scenario domains of char- 
acteristic size are formed, each of them has a random 2 . 5 
orientation chosen with uniform probabihty distribution 
from the D dimensional sphere. With a bias field this 



probability is not uniform. The landscape just after a 
transition is a uniform sea of spin down with islets of 
spin up scattered here and there. With an appropriate 
winding of spin in its neighbourhood an islet can survive 
as a localized topological texture. -7 
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Fig.l. The phase x (compare Eq.(^) as a function 
of X for an instantaneous quench with B ~ 0.005 and 
T = 0.00035. The thin line is the phase at t = 10 and 
the bold line is a plot for t = 1000. A (baby) Bloch wave 
can be identified as a place where one of the —2tt, 0, +2tt 
gridlines is crossed. 
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Fig. 2. Average density of Bloch waves 

as a function of temperature for magnetic fields: 
0.00125,0.0025,0.005,0.01 - from left to right. Each 
point is an average over several tausend runs with rela- 
tive error of 0.3% (less than the points size). The lattice 
constant was 0.2 and the periodic lattice size was 1000, 
which is 30 . . . 100 times the second correlation length. In 
each run the solitons were counted at i = 5/B (5 times 
the second relaxation time), when the soliton pattern is 
already stable - compare the bold line in Fig.l. Nimrer- 
ical data for B ~ 0.005 were fitted with the curve (HJ 



(solid line) to yield coefficients Ci = 0.88 — 0.02 and 

+ 

C2 = 2.45 - 0.04. The errors result from the fit. Other 
solid lines are extrapolation in good agreement with the 
data. 
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Fig. 3. Density of Bloch waves as a function of the 
quench time r for the fixed magnetic field B — 0.005 and 
the temperature T — 0.00035. The sohtons were counted 
at t = 2t + 5/B. Numerical data were fitted with the 

curve ( p^ - solid line. The coefficients are Di — 0.36 — 

0.005 and D2 = 0.74 - 0.02. 
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